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ABSTRACT 


The  pressure  and  phase  distribution  of  sound  in  a  fast  fluid 
medium  underlying  a  tapered  fluid  medium  was  modeled  using  a  Green's 
function  approach.  The  model  predicted  well  defined  beams  in  the 
bottom,  with  apparent  interference  effects  strongly  evident  for  cases 
at  a  distance  from  the  wedge  interface.  The  effects  of  attenuation  on 
the  patterns  were  studied.  A  simple  expression  for  determining  the 
beam  depression  angle  as  a  function  the  wedge  angle,  the  sound  speed 
ratio,  oUid  the  density  ratio  was  derived.  Comparisons  with  laboratory 
measurements  were  made. 
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I .  INTROJUCTION 


The  problem  of  radiation  behavior  in  a  wedge-shaped  medium 
overlying  a  bottom  possesing  a  greater  speed  of  sound  has  been 
investigated  both  optically  and  acoustically.  In  1971,  Tien  and 
-Martin  (Ref.  1)  deposited  a  thin,  tapered,  dielectric  film  on  a 
substrate  having  a  higher  refractive  index,  and  observed  the  behavior 
of  a  laser  beam  coupled  into  the  film.  Perfect  reflection  at  the 
film-substrate  interface  was  observed  until  the  changing  angle  of 
incidence  decreased  below  the  critical  angle  for  continued  mode 
propagation.  Por  angles  of  incidence  below  critical,  the  light  was 
converted  into  radiation  modes  in  the  substrate. 

A  concurrent  analysis  of  the  acoustic  problem  was  done  by  Kuznetzov 
(Ref.  2).  Kuznetzov  developed  a  theory,  based  on  normal -mode s ,  that 
predicted  the  following  for  sound  propagating  within  a  wedge-shaped 
medium  overlyi.ng  a  half-space  with  a  higher  sound  speed: 

1.  As  the  sound  approached  the  vertex  of  the  wedge,  it  would 
continue  to  be  totally  reflected  until  the  angle  of  incidence  decreased 
to  the  limiting  angle  of  total  reflection. 

2.  Sound  incident  at  less  than  the  limiting  angle  of  total 
reflection  would  be  totally  refracted  into  the  underlying  half-space. 
This  refraction  occurs  within  the  region  from  the  wedge  apex  to  the 
point  where  the  limiting  angle  was  attained. 


3.  Within  the  half -space ,  the  acoustic  energy  would  he  propagated 
as  a  well  defined  beam  with  the  maximum  pressure  amplitude  occurring  at 
an  angle  of  depression  from  the  wedge  interface  within  the  range  i  to  2B, 
where  e  is  the  wedge  angle. 

Kuznetzov  also  performed  a  series  of  experiments  that  supported  his 
theory . 

Subsequent  optical  experiments  by  Tien,  Smolinsky ,  and  karti.n  (Ref.  3) 
demonstrated  the  formation  of  a  well-defined  beam  within  the  substrate . 

The  accompanying  theoretical  discussion  predicted  conflicting  behavior . 
Ray-optics  predicted  refraction  into  the  substrate  after  the  cutoff  dis¬ 
tance  (the  point  where  refraction  begins)  and  wave  theory  predicted  com¬ 
plete  refraction  into  the  substrate  before  the  cutoff  point.  Sigelman , 
et.  al,  (Ref.  4)  both  predicted  a.nd  observed  a  well  columnated  beam  in  a 
water-aluminum  wedge  system,  but  the  observed  beam  was  much  broader  than 
predicted.  Also,  the  pressure  maxima  for  the  columnated  beam  occurred  at 
approximately  11°  far  a  wedge  angle  of  1.3°,  well  outside  the  limit  pre¬ 
dicted  be  Kuznetzov . 

The  primary  objective  of  this  theoretical  development  is  to  produce 
a  simple,  fast  computer  model  to  predict  the  pressure  and  phase  distri¬ 
bution  of  the  columnated  sound  beam  in  the  fast  substrate.  The  pressure 
and  phase  distribution  along  the  wedge  interface  calculated  by  Kawamura 
(Ref.  5)  will  be  used  as  input  to  the  substrate  model,  and  will  specifi¬ 
cally  be  used  as  the  boundary  conditions  for  a  Green's  function  analysis . 
To  remove  explicit  frequency  dependence ,  all  relevant  distances  will  be 
normalized  to  the  distance  measured  from  the  apex  along  the  wedge  inter¬ 
face  at  which  the  critical  reflection  angle  is  first  exceeded;  this  dump 


distance  is  referred  to  as  the  dump  distance  X.  Integration  over  the 


wedge  interface  will  he  truncated  at  distances  sufficient  to  encompass 
only  selected  modes  of  interest .  Results  will  be  analyzed  for  trends 
and  compared  with  e^ci sting  measured  data. 
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II.  THEORY 


Pressure  and  phase  distributions  along  the  bottom  of  a  wedge- 
shapad  medium  are  obtained  using  an  adaptation  of  the  method  of  Lmages 
as  prescribed  in  Ref.  5.  The  source  may  be  set  at  either  finite  or 
infinite  distance,  but  must  be  maintained  sufficiently  far  from  the 
wedge  boundaries  to  ensure  t.hat  impinging  waves  are  essentially  legally 
planar.  The  following  theoretical  discussion  uses  the  above  pressure 
and  phase  distribution  as  initial  conditions,  and  calculates  the 
pressure  and  phase  distribution  in  the  fast  bottom  by  utilising  the 
Green's  function  (Ref,  6). 

Figure  II-l  illustrates  the  geometry  of  the  problem.  The  horizontal 
x~a,'<is  and  the  vertical  t-axis  have  their  point  of  origin  at  the  wedge 
apex,  with  the  positive  x-axis  forming  the  interface  between  the  wedge 
and  fast  hottom.  The  wedge  angle  is  designated  3  and  the  angle 
specifying  the  sound  source  location  is  3.  X  is  the  first  mode 
dumping  distance,  and  is  a  distance  determined  as  follows-. 

3  arccos  (  C’  /  C’  )  (Il-la) 

c  ~  t  . 

0  =  Ct  /  {  4-f*sin(0  1  '  (Il-lb) 

0  ^  c 

X  =  D  /  sin(3)  (II-lc) 

0 
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where 


C[  =  sound  speed  in  the  wedge 

C;  =  sound  speed  in  the  bottom 

r  =  frequency 

and  =  lowest  mode  cut-off  depth. 

To  make  the  problem  frequency  independent,  all  distances  were 
normalized  to  X  by  defining 


-;  =  X  /  :<  (II-2a) 

=  X'  /  X  (Ii-2b) 

and  =  z  /  X  (i:-2c) 

The  variable  x'  is  identical  to  .x,  and  I'*  is  identical  to  .  The 
prune  designation  is  required  for  dummy  variables  of  integration  in 
the  Green's  function  development  to  follow. 

Pressure  along  the  wedge  bottom  is  then  given  by 

p(-;,r=0,t)  =  P(f)  exp  "iLwt  +  Ld)]'  (11-2) 


where  P(-;)  is  the  pressure  amplitude  along  the  wedge  bottom  and  l(') 
is  the  phase  along  the  wedge  bottom.  Pressure  in  the  bottom,  awa'/ 
from  the  boundary,  is  then  giver,  by  Ref,  7, 


p(r,t)  =  exp[iwt] 


P(-') 


X  dl' 


(II-4) 


z'=0 


=  exp  [i'JJt]  R(r) 


where  G  is  the  appropriate  Green's  function.  The  z'  coordinate  is 
identical  to  the  z  coordinate  and  r  =  (x-  +  z^)  ' .  The  appropriate 
Green's  function,  for  a  line  source  in  cylindrical  geometry  is 

G(^,r')  =  iTH^-^  (k2:?  -r''}  (II-5) 

where  k2  is  the  propagation  vector  in  the  bottom  and  r  =  [ (x')“  +  iz') 
Differentiation  of  the  Green's  function  yields  (P.ef.  7) 

^  =  iTH<"-'  (ko|r  -  -  (II-6) 

3z  '  - '  ■  i-'  t 

r  -  r  ' 

where  ’<2 1 1  "  r' ;  >>  Gtt  and  ,  H'^are  zero  and  first  order  Bessel 
functions  of  the  third  kind  (Hankel  functions)  respectively.  Using 
equations  (II-2a,b,c)  at  z''  =  0  the  expression  ir  -  r' '  may  be 
rewritten  as 


(II-7) 


Substituting  (II-7)  into  equation  (II-6) ,  and  evaluating  at  z'  =  0 
yields 


Combining  equations  (II-8)  and  (11-4)  yields  a  time  independent 


expression  for  amplitude  and  phase,  R(r) 


Equation  (II-9)  is  the  basic  equation  used  for  the  calculation  of 
phase  and  amplitude  in  the  bottom. 
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III.  DEVELOPMENT  OP  WDGBTM 


Program  WDGBTM  was  designed  to  perform  four  main  functions. 

First,  the  pressure  and  phase  along  the  interface  of  a  wedge  and 
underlying  fast  media  may  be  calculated  in  two  ways,  for  both 
infinitely  and  finitely  distant  sound  sources  in  accordance  with 
Ref.  5.  Modifications  were  made  to  the  code  of  Ref.  5  to  facilitate 
subsequent  in-bottom  calculations,  and  to  enhance  user  flexibility. 

These  changes  did  not  affect  results .  Then  using  the  first  code  as 
input,  additional  code  enaJales  both  the  computation  of  pressure  and 
phase  along  a  line  perpendicular  to  the  wedge  bottom,  and  as  a  function 
of  radius  and  angle  measured  from  the  first  mode  dumping  point,  X. 

Modes  of  program  execution  may  be  chosen  by  literal  input  directives 
combined  with  appropriate  numerical  input.  Appendix  A  illustrates  an 
input  set  that  utilizes  all  four  options.  Output  consists  of  both 
printed  listings  and  graphics. 

Program  versions  exist  for  execution  on  both  a  CDC  6500  and  an  IBM 
360.  The  CDC  6500,  with  a  60  bit  word,  does  not  require  double 
precision  for  accuracy  equivalent  to  the  IBM  360  output.  The  CDC 
version  of  WDGBTM  has  simpler  varian  control  routines,  and  separate 
output  files  for  each  pressure/phase  variation.  This  enables  optional 
suppression  of  listings.  Execution  time  on  the  CDC  6500  is  roughly  ten 
times  faster  than  on  the  IBM  360.  For  average  input  sets,  times  are  150 
and  1500  cp  seconds  respectively. 
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A.  EQUATION  PREPARATION 
1 .  Depth  Domain 

Integrations  along  lines  perpendicular  to  the  wedge  bottom 
were  designated  depth  domain  calculations.  The  program  was  designed  to 
allow  specification  of  the  field  point,  i.e.,  specification  of  where 
along  the  media  interface  line  the  perpendicular  probe  line  could  be 
dropped.  The  field  point  is  not  confined  to  lines  directly  beneath  the 
wedge,  but  may  extend  to  points  beyond  the  wedge  apex.  These  points 
would  be  values  of  negative  '■<  in  Figure  ll-l.  The  beginning  and  ending 
points  of  the  section  of  the  perpendicular  to  be  examined,  and  the 
number  of  field  points  used  in  the  interval  were  also  to  be  specified. 

Equation  (11-9)  was  split  into  real  and  imaginary  components  before 
coding.  Integration  order  selected  was  over  the  entire  5 'domain  for 
each  n  point. 

2 .  Angle  Domain 

Examination  of  Figure  ll-l  verifies  the  following  geometric 
relationships . 


and 


£  =  0  COS0, 
k 


(iri-la) 


n  =  (o^  -  (iiz-ib; 

where  £  is  the  distance  from  X  to  the  5'  coordinate  of  the  field  point 
of  interest  and  0  is  the  radial  distance  from  X  to  the  field  point. 
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Substitution  of  equations  (Ill-la)  and  (Ill-lb)  into  equation  (II-9) 


yields 


/  ^  -2,1/2 

R(r)  =  — — — -  exp 


/  in  \  I2k-,X\ 


1/2 


P  (Z ')  exp{i[L(Z ') 


>0 


(III-2j 


Equation  (111-2)  is  the  basis  for  angle  domain  calculations.  The 
specification  of  a  originating  at  X  was  chosen  to  conform  to  e.xisting 
experimental  data  (Reference  3).  The  length  of  the  radius,  the 
beginning  and  ending  of  the  angular  sweep,  and  the  angular  increment 
comprise  the  input  specifications.  As  in  the  depth  domain  option, 
equation  (III-2)  was  split  into  real  and  imaginary  components  before 
coding . 


3 .  Design  Restrictions 

Theoretical  development  of  program  WDGBTM  assumes  the  sound 
source  is  sufficiently  distant  to  ensure  essentially  locally  planar 
waves  at  the  interface  and  the  in-bottom  measurements  satisfy  the 
condition  |t-rJ>>2'rT. 

4 .  Program  Results 

An  initial  set  of  WDGBTM  runs  were  made  in  the  depth  domain  for 
a  perpendicular  directly  dov/n  from  the  wedge  apex.  The  source  was 
infinitely  distant  and  set  at  an  angle  of  B/2.  This  limited  modes  within 
the  wedge  to  odd  numbered  (symmetric)  inodes  only.  Integration  over  the 
wedge  bottom  was  truncated  at  %  equal  four,  so  that  only  modes  one  and 
three  were  considered  as  contributing  significantly  to  the  beam  pattern 
in  the  bottom.  Beam  patterns  in  the  bottom  were  calculated  over  the  range 
h  equal  zero  to  1.3X. 
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The  initial  set  was  comprised  of  seventy- five  runs.  In  addition 
to  the  above  criteria ,  the  runs  were  generated  by  all  possible  combi¬ 
nations  of  the  following  specifications.  The  ration  ^ 2  0.95, 

0.90,  0.80,  0.70  or  0.50.  The  sound  speed  ratio,  C^/C^  was  set  at  0.95, 
0.90,  or  0.85.  The  wedge  angle  3  was  varied  through  10,  6.92,  5,  2  and 
1  degree (s)  . 

From  the  results  of  the  initial  data  set,  the  depression  angle  and 
beam  width  was  calculated  for  t.he  beam  pattern  in  the  bottom.  Results 
are  summarized  in  Appendix  B.  These  results  were  then  used  to  find  a 
simple  analytical  formula  to  predict  the  depression  angle  as  discussed 
in  Section  IV. 

Other  runs  were  made  for  perpendiculars  at  x  =~9Xfrom  the  wedge 
apex.  These  runs  displayed  apparent  strong  interference  phenomena  and 
are  discussed  in  detail  later. 

Angle  domain  executions  were  made  to  compare  with  experimental  data. 
Results  are  discussed  in  Section  VII. 

a.  SENSITIVITY  TO  INTEGRATION  STEP  SIZE 

It  was  desirable  to  keep  both  the  integration  step  size  and  the  field 
point  increment  size  as  large  as  possible  in  order  to  minimize  total 
computer  resources  required .  Increment  of  0.01  times  the  dump  distance 
(X)  proved  to  be  optimum  for  both  integration  increments  along  the 
wedge  interface  and  field  point  steps  in  the  bottom. 

To  illustrate  the  adequacy  of  O.OIX  increment,  a  series  of  runs  were 
made  in  accordance  with  Table  III-l.  An  overlay  of  wedge  interface 
pressure  amplitude  for  lx  =  O.OIX  and  lx  =  O.OOlX  demonstrated  no  dis¬ 
tinguishable  difference.  Likewise,  depth  domain  plots  (a,  b,  c,  d  of 
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Table  (III  -  1) 

Sensitivity  to  Integration  Step  Size  (i-v)  and  Field  point  interval 


size  (An) 


A-v  =  0.001 


An=o.(:i  An=o.ooj 


n  amp  max 


0.240 


0.244 


0.240 


0.243 


iY  =  0.0001 


■lr\=0.01  Art=0.001 


13.71 


13.50 


13.66 


Table  could  be  identically  superimposed .  Use  of  an  integration 

step  size  of  O.OOIX  resulted  in  a  calculated  depression  angle  within 
one  percent  of  the  depression  angle  for  the  O.OlX  run.  Utilization  of 
a  step  size  of  greater  than  O.OlX  is  not  recommended.  Within  this 
limitation  the  integrations  can  be  carried  out  over  lOX  along  the  inter¬ 
face  for  a  maximum  field  interval  of  lOX  in  the  bottom.  The  lOX  re¬ 
striction  is  imposed  by  program  array  size.  (For  the  T‘s  specified ,  0.5 
degrees  in  the  angle  domain  was  slightly  finer  than  the  O.OlX  field 
point  step  size  in  the  depth  domain  and  was  the  field  point  increment  of 
choice) . 
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IV.  BEAM  DEPRESSION  ANGLE  ANALYSIS 


■A  simple  relationship  was  sought  between  the  physical  parameters 
associated  with  the  wedge  system  and  the  resultant  beam. 

Assume 


-DC 


(IV-1) 


where  -i  is  the  beam  depression  angle  and  c  is  the  wedqe  angle.  For 
DC 

each  and  Cj^/C,  ,  ®  graphed  on  log-log  paper,  and  the 

slope  {  n)  was  measured.  Results  are  listed  in  Table  1  of  Appendix  C, 
n  equals  0.329. 

Then  >5^ was  found  as  a  function  of  The  critical  angle,  , 

is  an  explicit  function  of 


=  arccos  (c,/C_)  in  radians 

L  12^ 


It  was  then  assumed 


(IV-2) 


,0. J29 


DC 


(IV-3/ 


For  any  constant  3 ,  .V  was  derived  as  the  slope  of  a  log-log  graph  of 
vs  0^.  Graphs  were  analysed  for  3  equal  to  10,  6.92,  5,  2  and  1 
degree,  and  ratios  of  0.95,  0.90,  0.35,  0.70  and  0.50.  Results 

are  summarired  in  Table  (2)  of  Appendix  C,  ^  =  0.772. 
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The  same  process  was  repeated  for  ascertaining  dependence  on 


.  Log-log  graphs  were  prepared  by  holding  3  and  constant 
using  the  assumed  relation 


^  ^  .0.329  ,  0.772  ,,  ,  ,g 


(IV-4) 


Results  are  summarized  in  Table  3  of  Appendix  C.  q  =  -0.254. 

The  proposed  formula  is  now  dependent  on  all  physical  properties, 
assuming  the  source  angle  is  maintained  at  one-half  3.  To  remove  the 
proportionality  sign,  write 


.0.329  ,  0.772 


('0 


-0.254 


( IV-5 ) 


and  solve  for  K  .  From  Table  4,  Appendix  C,  =  17.22/radian. 
Therefore 


,  ^0.329  .  0.772  -0.254 

0^  =  17.22  3  0,  (0^/0^) 


(TV-6) 


where  3  is  in  degrees  and  ^(-in  radians. 

a„  calculated  from  equation  (IV-6)  and  the  difference  between  0  _ _ 
CL  CL 

and  0^  from  program  WDGBTM  (labeled  and  0  respectively )  are  also 
tabulated  on  Tcible  4,  Appendix  0.  (10^  varies  from  a  minimum  of  0.00  to 

a  maximum  of  1.25,  with  the  average  of  the  absolute  values  of  equal 
to  0.36.  Errors  ranged  from  zero  percent  to  5.0  percent.  Thus  gives 
a  reasonable  approximation  of  the  beam  depression  angle. 


A  similar  formula 


iIV-7) 


was  developed  independently  by  A.  Coppens  (Reference  6)  The  minimum 

deviation  of  d _ from  9„.is  0.01  and  the  maximum  deviation  is  1.56. 

DA  DM 

The  average  of  the  absolute  value  of  the  deviation  is  0.36.  Both 
formulas  are  best  for  P-^/O^  close  to  1.0,  and  less  accurate  for 
0-^/p2  ®<iual  to  0.5. 
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V.  THE  EFFECT  OF  ATTENUATION 

As  previewed  in  Section  III~A.4  the  pressure  distribution  as  a 
function  of  0  at  -v  =  -9. OX  demonstrated  apparent  interference  effects 
(see  figure  V-1) .  As  integration  over  the  wedge  interface  to  4X 
encompasses  dumping  by  the  first  and  third  mode,  and  the  third  mode 
dumps  about  the  point  3X  vice  X  for  the  first  mode,  it  is  log'  i.  to 
expect  attenuation  effects  to  reduce  the  contribution  of  the  more 
distant  mode.  The  initial  Green's  function  approach,  equation  (II-9) 
did  not  consider  attenuation  of  the  acoustic  wave  in  the  bottom 
medium.  An  attenuation  term 


exp 


(V-l) 


(V-:) 


’  '  -J' 


was  incorporated  into  equation  (II-9)  to  yield 
R(r)  =  n/4  exp^^J 

';exp''-a  ‘-r- j  ‘^‘■'e.Yp{i  rii'i;  ■ 

An  appropriate  code  modification  to  the  depth  domain  routine  was  made. 
The  attenuation  coefficient  ->•  was  an  input  variable  to  the  program. 

Two  sets  of  comparison  runs  were  made.  The  first  set  consisted  of 
an  unattenuated  control  (9  =  0)  and  an  attenuated  run  with  -  equal  1.0, 
with  the  pressure  amplitude  measured  down  from  the  wedge  apex.  The 
second  set  consisted  of  an  unattenuated  control  plus  runs  with  ^  set  to 
1.0,  0.5,  0.1,  and  0.01  and  the  pressure  measured  at  =  -9.0  from  the 
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apex.  All  other  input  variables  for  both  sets  were  held  constant  and 
are  specified  in  Table  (V-1) .  Integration  over  the  wedge  bottom  was 
from  the  apex  to  four  times  the  dump  distance;  this  was  to  insure  that 
contributions  from  the  third  mode  (dumping  about  the  point  x=  3X)  were 
included. 

The  results  of  set  one  are  illustrated  in  Figure  (V-2) .  Note  that 
the  second  maxima  is  severely  depressed  in  the  attenuated  case.  The 
depression  angle  has  been  shifted  slightly  from  11. S6  degrees  for  the 
unattenuated  data  to  11.31  degrees  with  x  set  to  1.0. 

Figure  V-3  is  a  plot  of  the  results  from  set  two.  All  curves  have 
been  normalized  to  the  initial  slope  for  =  1.0.  The  second  and  third 
maxima  decrease  with  increasing  d,  and  relatively  more  energy  is 
contained  in  the  first  beam.  The  depression  angle  of  the  first  beam 
decreases  slightly  with  increasing  a  .  These  trends  are  summarized  in 
Table  {V-2). 

Figures  (V-4a)  and  (V-4b)  illustrate  the  effects  of  attenuation  on 
different  contributing  modes.  For  Figure  (V-4a) ,  the  source  was  placed 
at  1/3  S.  With  this  source  geometry,  contributions  from  mode  three  are 
suppressed,  and  only  the  first  and  second  mode  contributed  to  the  beam 
pattern  in  the  bottom.  For  Figure  (V-4b) ,  the  source  was  placed  at  43 
and  modes  one  and  three  propagate  into  the  bottom  medium.  Both  cases 
were  run  for  u=  1.0,  0.5,  0.1  and  0.0.  As  could  be  expected, 
attenuation  causes  greater  suppression  of  third  mode  interference 
contributions  than  of  second  mode  contributions. 
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Table  (V-1) 


Input  Specif ications  for  Attenuation  Runs 


o 

tl 

1500.0 

=  1.0 

'^2 

1578.947 

=  1.052632 

2 

Ci/c^ 

=  0.95 

P^/O^  =  0.95 

S 

O 

O 

in 

It 

f  =  150000  hz 

^2 

6 

o 

o  m 

II  II 

Source  distance  = 

a 

tl 

o 

o 

o 

o 

0.1,  0.5,  1.0 

Refer  to  Figure  II-l  for  variable  meaning. 
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Figure  V-2 


Attenuated  and  unattenuated  pressure  amplitudes 
down  from  the  apex 


o.tf  ase  e.rs  i.oo  i.tr 


ETft  - 


26 


TABLE  V-2 


1 

Y 

H^/H 

1.0 

12.84 

0.21xl0~‘^ 

NONE 

0.5 

13.39 

.33xlO~^ 

.26xl0~^ 

1.26 

0.1 

13.82 

0.21 

0.23 

0.90 

0.01 

13.93 

0.54 

0.66 

0.  83 

0.0 

13.98 

0.61 

0.74 

0.82 

Y  =  depression  angle  of  first  maxima 
=  amplitude  of  first  maxima 
=  amplitude  of  second  maxima 
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VI.  ANALYSES  OF  FAR-FIELD  INTERFERENCE  PATTERNS 


The  results  of  the  attenuation  investigation  were  supportive  of  an 
interference  explanation  for  the  multiple  beams.  Subsequently,  two 
analytical  approaches  were  tried  to  predict  the  location  of  interference 
maxima  and  minima  under  unattenuated  far-field  conditions. 

Both  the  Differential  Change  and  Multiple  Lambda  analyses  assume 
very  simple  processes.  Both  assume  point  sources  of  sound  at  the  wedge 
interface  loc^ited  at  X  and  iX.  Reallistically  the  sound  sources  behave 
as  phased  arrays  propagating  assymetrically  about  the  dumping  points. 
Phase  differences  between  the  two  dumping  points  are  non-zero.  The 
Differential  Change  formula  becomes  invalid  as  n  approaches  zero. 

A.  DIFFERENTIAL  CHANGE  APPROACH 

Based  on  the  simple  geometric  considerations  illustrated  in  Fig  (VI-1) 
the  following  nondimens ional  formula  for  the  change  in  interference 
maxima  locations  was  derived. 

4  (C^/C^)  sin  0^  sin  S 

^  ~  .  "f  (VI-1) 

Note  that  equation  (VI-1)  is  dependent  only  on  the  wedge  angle  3  and  the 
sound  speeds  of  the  two  fluid  media. 

A  series  of  WDGBTM  executions  were  made  for  varying  C^/C^  and  with 
all  other  variables  held  constant.  Table  (VI-1)  is  a  summary  of  one  case 
this  case  is  typical  of  results  observed.  Results  will  be  discussed  in 


TABLE  VI-1 


Comparison  of  An  from  equation  (VI-1)  and  An  from  WDGBTM  results. 


i  WDGBTM 

maxima 
locations 

3.0  2.38 

4.05 
5.  97 

2.5  2.33 

3.73 
4.75 
5.66 

1.25  2.28 

3.30 
3.96 
4.45 


C  /C^  =  0.90 


-in 

from 

WDGBT.'^ 


An 

predicted 


1.67 
1 . 92 


4 . 56 
3.04 


1.40 

1.02 

0.91 


2.  33 
1.60 
1.38 


1.22 
0.  66 
0.49 


1.29 
0.88 
0.  77 
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B.  MULTIPLE  LAMBDA  APPROACH 

The  second  interference  pattern  calculation  also  utilizes  the 
geometry  of  Fig.  (Vl-1) .  For  each  ratio  and  beta  combination, 

X  and  (the  wavelength  in  the  bottom  media)  was  calculated  according 
to  the  equations 


(VI-2) 

(VI-3) 

where  is  ISOkhz.  The  ratio  X^/x  =  Q 
is  then  found  and  used  in  the  relation 


X  = 


4f  sin  9^  sin  6 


and 


X  =  c  /f 
2  2' 


w  =  (vi-4) 

where  X  is  any  positive  integer.  Note  that  use  of  equation  (VI-4) 
assumes  zero  initial  phase. 

To  find  the  expected  maxima  and  minima,  a  plot  of  vs  ^  was  made 
for  each  Q.  Maxima  are  expected  for  each  integer  N,  minima  for  each 
.V  +1/2.  Tcible  (VI-2)  summarizes  one  case  of  maxima  and  minima  expected 
as  a  result  of  equation  (’1-4)  versus  maxima  and  minima  observed  on 
WDGBTM  output. 

C.  SUMMARY  OF  ANALYSES 

Predictions  from  both  simple  analyses  are  of  the  Scime  order  of 
magnitude  as  WDGBTM  results.  The  number  of  maxima  and  minima  are 
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TABLE  VI- 2 


OCCURBEtlCE  OF  MAXIMA  AND  MINIMA 
C^/C^  =  0.3S 


2 {DEG) 

n  MAX  WDGBTM 

n  MAX  EXP 

n  min  mdgbtm 

n  MIN  EXP 

5 

4.03 

2.68 

5.  39 

4.70 

6.34 

2.  5 

3.  59 

2.59 

4.  55 

3.  76 

5.00 

4.70 

5.  66 

5.56 

2.25 

2.02 

2.53 

2.49 

1.76 

2.  35 

3.  56 

3.16 

3.31 

3.72 

4.77 

4.25 

5.03 

4.68 

5.  40 

5.13 

5.62 

5.54 

5.  90 

5.  94 
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consistantly  predicted  as  increasing  with  decreasing  wedge  angle. 
Occasinal  unexplainable  phenomena,  such  as  increasing  WDGBTM  In  for 
the  3  =  5  on  Table  (VI-1)  and  the  extra  WDGBTM  maxima  f  :)r  3=5 
Table  (VI-2),  were  observed,  but  results  are  generally  supportive  of 
an  interference  effect  occurring. 
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VII.  COMPARISONS 


A.  EXPERIMENTAL  DATA 

Netzorg  (Reference  4)  measured  pressure  amplitude  as  a  function  of 
angle  for  fixed  radii  in  four  media  below  a  wedge  shaped  upper  layer. 

The  angle  domain  version  of  WDGBTM  was  used  to  compare  calculated 
pressure  patterns  with  measured  patterns. 

All  runs  utilized  finite  source  input.  For  each  of  Netzorg 's  four 
experimental  cases,  corresponding  computer  runs  were  made  with  the 
source  distance  set  to  lOOX,  lOX  and  SX,  where  S  corresponds  to  the 
actual  source  distance  used  in  the  tank.  In  all  cases  S  was  less  than 
lOX.  The  integrations  were  carried  out  over  one  dumping  distance  and 
four  dumping  distances  for  each  source. 

Figure  (VII-1)  is  typical  of  pressure  amplitudes  in  the  bottom. 

The  source  in  the  experiment  was  only  about  40  wavelengths  from  the 
dumping  point,  even  at  lOX  the  source  is  approximately  IIC  wavelengths 
away.  Since  the  method  of  images  calculations  utilized  to  determine 
the  pressure  aunplitude  along  the  wedge  bottom  assume  a  plane  wave,  it 
is  questionable  that  this  program  can  produce  a  realistic  comparison 
with  the  measured  data. 

Figures  (VII-2)  and  (VII-3)  illustrate  the  wide  variation  in  results 
obtained  depending  on  the  source  distance  and  integration  interval. 
Calculations  demonstrate  expected  trends,  i.e.,  pressure  amplitude 
increases  with  decreasing  radius  distance  (") ,  with  closing  source 
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distance,  and  as  a  function  of  an  increasing  integration  interval  over 
the  wedge  bottom.  The  beam  depression  angle  becomes  more  shallow  as 
source  distance  increases.  Results  at  very  low  angles  are  not 
considered  valid  because  of  assumptions  inherent  in  the  WDGBTM  model. 

The  circled  data  points  of  Figures  (VII-2)  and  (VII-3)  are  a 
transposition  of  data  from  Reference  4.  Since  absolute  amplitude  was 
not  specified  in  the  reference,  rescaling  to  the  case  most  near  to  the 
actual  conditions  was  performed,  i.e. ,  the  data  were  normalized  to  the 
case  =  1.0  and  srcd  =  3.5X.  The  shape  of  the  experi.mental  curve  and 
the  apparent  beam  depression  angle  are  the  only  valid  comparisons  that 
should  be  made. 

Comparisons  of  WDGBTM  beam  depression  angles  with  measured  values 
gave  varying  results.  Measured  depression  angles  decreased  with 
increasing  J  ,  WDGBTM  depression  angles  increased.  WDGDTM  results  were 
within  6  percent  to  62  percent  of  measured  values.  No  correlation  was 
observed  between  prediction  accuracy  and  source  distance  or  prediction 
accuracy  and  the  interface  integration  interval.  All  measured  beam 
depression  angles  in  Reference  3  were  for  the  two  shortest  radii  of 
each  case,  the  three  longer  radii  could  only  be  qualitatively  compared 
with  graphical  data.  Possible  attenuation  effects  were  not  considered. 
Numerical  results  are  contained  in  Appendix  D. 

B.  PARABOLIC  EQUATION  APPROACH 

Kuperman  (Reference  9)  describes  a  Parabolic  Equation  model  result 
for  a  wedge  problem,  and  compares  these  results  to  one  set  of  data 
extracted  from  Reference  10. 
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The  Parabolic  Equation  results  show  very  well  defined  beaims  in  the 


bottom,  with  no  apparent  interference  effects  as  observed  in  comparable 
WDGBTM  runs  (Figures  V-1  and  V-4) .  Kuperman  used  an  attenuation  factor 
of  0.5dB  per  wave  length,  a  value  equivalent  to  an  input  of  0.984 
Nepers/X  for  WDGBTM.  It  should  be  noted  in  Figure  V-4  that  there  is 
complete  separation  of  the  beams  for  WDGBTM  when  ^  =  1.0. 

Kuperman  states  agreement  within  20  percent  of  a  measured  data  set 
of  Reference  10.  This  measured  data  is  identical  to  Case  1  of  Appendix 
D.  WDGBTM  accuracy  for  this  case  varied  from  zero  percent  error  to  22 
percent  error  with  error  calculated  according  to  the  formula 

^  (observed  data)  -  (measured  data) 
measured  data 
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VIII.  CONCLUSIONS 


1.  Progrcun  WDGBTM  offers  a  fast,  efficient  means  of  analysing  the 
behavior  of  sound  in  a  fast  bottom  underlying  a  wedge-shaped  medium. 

Effects  of  attenuation  are  easily  demonstrated,  and  modes  desired  are 
readily  controlled  by  specification  of  source  location  and  integration 
interval  selected.  Results  of  the  model  are  not  inconsistant  with 
current  measured  data,  and  illuminate  observations  generated  by  the 
Parabolic  Equation  model. 

2.  A  single  equation  for  calculating  beam  depression  angle  as  a 
function  of  wedge  characteristics  was  obtained  which  offers  sufficient 
accuracy  for  use  as  a  rapid  analytical  tool.  The  equation  was  derived 

for  acoustic  sources  placed  at  one-half  the  wedge  angle.  Accuracy  degrades 
when  the  wedge  medium  density  to  bottom  medium  density  ratio  decreases 
to  0.5. 

3.  It  is  recommended  that  further  experiments  be  performed  that 
would  further  test  the  model.  New  experiments  should  be  sufficiently 

large  scale  to  ensure  plane  waves  at  the  wedge  bottom  interface  and  to  satis- 
fy  the  k^\r-r  criteria.  Bottom  media  that  reasonably  model 

real  ocean  fluid-like  bottom  are  also  recommended. 


APPENDIX  A 


SAMPLE  INPUT  SET 


o 

o 

O 

o 

o 

<3 

5.0 

>o 

'O 

>0 

• 

OsJ 

in 

C<\ 

CO 

sO 

CO 

n 

00 

ou> 

o<^ 

•  « 

•  • 

»-^co 

m4««4 

o 

o 

« 

• 

o 

O' 

(T 

• 

» 

1 

o 

o  • 

OO 

>00 

OO 

>oo 

OO 

QOO 

OO 

OO 

m  in 

1) 

•m 

•v> 

•  • 

0-1 

o  o 

a 

c 

•H 

r*- 

0 

ro 

o 

o 

0 

0 

'-H 

o 

o 

•  *0 

— f 
r\ 

• 

• 

m  in  • 

U.4 

>c 

o 

^  • 

• 

c 

aoo 

Koo 

0 

zc^o 

►- 

Ocvio 

»-  >-  1- 

•-^tn  • 

o 

o 

a>o  • 

O  o  o 

Qj 

a 

a 

Z-->t 

a  0.  a 

UJ 

z 

z 

z  z  z 

u 

M  M 

0^ 

UJ 

o  o 

o 

z 

z 

o 

z  •z  •z 

a 

0^ 

N^0*-^0“ 

to 

< 

o 

<  <  < 

z 

z 

o 

z  z  z 

Ui 

o 

o 

VO 

o  o  o 

o 

o 

O  Q  O 

•-•o 

UjO 

o 

Z  t 

X 

ox 

OH-  • 

•  UJ  UJ  X  o 

— o 

OH- 

00-J-*-J0»-  o 

u.oooa  0-— a  oozooo-jO'^ooq.  oo 

•  *UJ  •  ••-irtlA  3  *Z  •LU  •  • 

-^•^O  QOOQOOU.-<0  <--<OQOO 


44 


Litoral  string  for  infinite  source  option 


APPENDIX  B 


3^/p2  =  -9500 


Cj/C^ 

B (DEG) 

“^DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

15.64 

22.93 

0.68 

6.92 

13.50 

19.08 

0.71 

5.00 

12.02 

17.48 

0.69 

2.00 

9.09 

12.19 

0.74 

1.00 

7.58 

9.54 

0.79 

0.90 

10.00 

19.80 

29.07 

0.68 

6.92 

17.74 

26.66 

0.67 

5.00 

16.33 

21.35 

0.75 

2.00 

11.31 

16.20 

0.70 

1.00 

9.09 

12.73 

0.71 

0.85 

10.00 

24.37 

36.02 

0.68 

6.92 

21.16 

31.05 

0.68 

5.00 

18.42 

26.61 

0.69 

2.00 

13.50 

19.45 

0.69 

1.00 

11.31 

14.36 

0.79 

-^DEPRESSION  =  ^CTAN(Z/X)  AT  MAX  AMPLITUDE 

BANDWIDTH  =  AECTAN(Z/X,  VX,  ) 

lU  -  ‘jL 

l£ 

MEAM  BW  =  .73 
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9000 


C^/Cj 

6  (DEG) 

4. 

^DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

15.64 

22.68 

0.69 

6.92 

13.49 

19.08 

0.71 

5.00 

12.02 

17.43 

0.69 

2.00 

9.09 

12.19 

0.74 

1.00 

7.58 

9.54 

0.79 

0.90 

10.00 

19.80 

29.77 

0.66 

6.92 

18.42 

26.75 

0.69 

5.00 

16.33 

21.80 

0.75 

2.00 

12.02 

15.53 

0.77 

1.00 

9.09 

12.90 

0.70 

0.85 

10.00 

24.37 

36.28 

0.67 

6.92 

21.16 

31.22 

0.68 

5.00 

18.42 

26.93 

0.68 

2.00 

13.50 

19.70 

0.68 

1.00 

11.31 

14.95 

0.76 
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8000 


C  0  =  . 

^ir2 


C1/C2 

S(DEG) 

■^DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

16.33 

22.39 

.73 

6.92 

14.20 

19.19 

.74 

5.00 

12.79 

17.33 

.74 

2.00 

9.82 

12.13 

.81 

1.00 

7.58 

9.87 

.77 

0.90 

10.00 

20.46 

31.00 

.66 

6.92 

19.14 

26.93 

.71 

5.00 

16.33 

21.90 

.74 

2.00 

12.02 

16.38 

.73 

1.00 

9.82 

12.95 

.76 

0.85 

10.00 

25.03 

36.94 

.68 

6.92 

22.44 

31.72 

.71 

5.00 

19.14 

27.65 

.69 

2.00 

14.20 

20.10 

.71 

1.00 

11.31 

15.32 

.74 
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7000 


C/Cj 

S  (DEG) 

4. 

DEPRESSION 

BANDWIDTH 

RATIO 

0.95 

10.00 

17.07 

22.24 

.77 

6.92 

14.20 

19.34 

.73 

5.00 

13.50 

17.17 

.79 

2.00 

9.82 

12.13 

.31 

1.00 

8.36 

9.65 

.87 

0.90 

10.00 

21.16 

32.66 

.65 

6.92 

19.80 

26.93 

.74 

5.00 

17.07 

22.00 

,78 

2.00 

12.79 

16,59 

.77 

1.00 

10.59 

13.12 

.81 

0,85 

10.00 

25.64 

37.88 

.68 

6.92 

23.12 

32.29 

.72 

5.00 

19.30 

28.72 

.69 

2.00 

14.95 

20.56 

.73 

1.00 

12.02 

15.85 

.76 
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5000 


r- 


^1/^2  ^ 


Ci/Cj 

3 (DEG) 

4- 

DEPRESSION 

BANDWIDTH 

K.ATXO 

0.95 

10.00 

17.74 

22.73 

0.78 

6.92 

15.64 

19.80 

0.79 

5.00 

14.95 

16.91 

0.88 

2.00 

9.82 

12.68 

0.77 

1.00 

9.09 

10.26 

0.89 

0.90 

10.00 

22.44 

35.53 

0.63 

6.92 

21.80 

27.16 

0.80 

5.00 

18.42 

23.70 

0.78 

2.00 

14.20 

17.64 

0.80 

1.00 

11.31 

13.87 

0.82 

0.85 

10.00 

27.20 

42.80 

0.64 

6.92 

25.03 

34.33 

0.73 

5.00 

21.80 

31.00 

0.70 

2.00 

17.07 

21.40 

0.78 

1.00 

13;  50 

16.91 

0.80 

I 


I 
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APPENDIX  C 


TABLE  1 


^^1/2 

Ci/C^ 

n 

0.95 

0.95 

0.318 

0.95 

0.90 

0.347 

0.95 

0.85 

0.341 

0.90 

0.95 

0.318 

0.90 

0.90 

0.335 

0.90 

0.85 

0.339 

0.80 

0.95 

0.328 

0.80 

0.90 

0.328 

0.80 

0.85 

0.343 

0.70 

0.95 

0.318 

0.70 

0.90 

0.303 

0.70 

0.85 

0.343 

0.50 

0.95 

0.339 

0.50 

0.90 

0.314 

0.50 

0.85 

0.314 

Average  n 

=  0.329 

Standard  (N-1) 

deviation 

=  0.014 

(N-1)  Variance 


0.002 


TABLE  2 


'^1/^2 

5 

0.95 

10.0 

0.834 

0.95 

6.92 

0.845 

0.95 

5.00 

0.783 

0.95 

2.00 

0.7  39 

0.95 

1.00 

0.752 

0.90 

10.0 

0.815 

0.90 

6.92 

0.818 

0.90 

5.00 

0.782 

0.90 

2.00 

0.735 

0.90 

1.00 

0.720 

0.80 

10.0 

0.791 

0.80 

6.92 

0.839 

0.80 

5.00 

0.767 

0.80 

2.00 

0.686 

0.30 

1.00 

0.750 

0.70 

10,0 

0.730 

0.70 

6.92 

0.892 

0.70 

5.00 

0,701 

0.70 

2.00 

0,778 

0.70 

1.00 

0.697 

0.50 

10.0 

0.786 

0.50 

6.92 

0.881 

0.50 

5.00 

0.693 

0.50 

2.00 

1.051* 

0.50 

1.00 

0.722 

Rejected  data 

point 

Average  Af 

=  0.772 

Standard  (N-1)  deviation 

=  0.059 

Variance 

=  0.003 
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Table  3 

Ci/c, 

Q^(rad) 

9 

0.95 

0.318 

10. 0 

'0.217 

0.95 

0.318 

6.^2 

'0.285 

0.95 

0.318 

5.00 

'0.257 

0.95 

0. 313 

2.00 

'0.215 

0.95 

0.318 

1.00 

-0.283 

0.90 

0.451 

10.0 

'0. 183 

0.90 

0.451 

6.92 

-0.232 

0 . 90 

0.451 

5.00 

-0.211 

0.90 

0-451 

2.00 

-0.300 

0.90 

0.451 

1.00 

-0.317 

0.35 

0.555 

10.0 

-0. 162 

0.85 

0.555 

6.92 

-0.242 

0.35 

0.555 

5.00 

-0.255 

0.95 

0.555 

2.00 

-0.362 

0.85 

0.555 

1.00 

-0.283 

Average  <J  =  -0.254 

standard  (N-1)  deviation  =  0.053 

variation  =  0.0028 


S3 


Table  4 

Ol/P  2 

C^/C2 

3 

^DW 

K 

®DC 

0 .95 

0 . 95 

10.0 

15.64 

17. 5u 

15.35 

-0 .29 

0.95 

0 .95 

6.92 

13.50 

17.08 

13 . 60 

0.10 

0.95 

0.95 

5.00 

12.02 

16,94 

12.22 

0 .20 

0.95 

0.95 

2 . 00 

9.09 

17,32 

9.04 

-0.05 

0.95 

0.95 

1.00 

7.58 

18.14 

7,20 

-0 .38 

0.95 

0.90 

10.0 

19,80 

16.94 

20.12 

0 . 32 

0.95 

0.90 

6.92 

17 . 74 

17.13 

17 .83 

0.09 

0.95 

0.90 

5.00 

16,33 

17.55 

16.02 

-0.31 

0.95 

0.90 

2 .00 

11.31 

16,43 

11.85 

0,54 

0.95 

0 .90 

1.00 

9.09 

16.59 

9.43 

0.32 

0.95 

0.35 

10.0 

24.37 

17.77 

23  .  61 

-0.76 

0.95 

0.85 

6.92 

21.16 

17. U2 

20.92 

-0 , 24 

0.95 

0.85 

5,00 

18,42 

16.87 

18.80 

0,38 

0.95 

0.85 

2.00 

13.50 

16.72 

13.91 

0.41 

0.95 

0.85 

1.00 

11.31 

17.59 

11.07 

-0.24 

0.90 

0.95 

10.0 

15.64 

17.31 

15.56 

-0.08 

0.90 

0.95 

6,92 

13.49 

16.85 

13.79 

0.30 

0.90 

0,95 

5,00 

12.02 

16.71 

12.39 

0.37 

0.90 

0.95 

2,00 

9.09 

17.03 

9.16 

0.07 

0.90 

0.95 

1.00 

7,58 

17.89 

7.30 

-0.28 

0.90 

0.90 

10.0 

19.80 

16.71 

20,40 

0.60 

0.90 

0.90 

6.92 

18.42 

17.55 

18.08 

-0 . 34 

0.90 

0.90 

5.00 

16.33 

17.31 

16 , 24 

-0.09 

0,90 

0.90 

2.00 

12.02 

17.23 

12.02 

0.00 

0.90 

0,90 

1.00 

9.09 

16.36 

9  .  56 

0.47 

0.90 

0.85 

10,0 

24.37 

17.53 

23.94 

-0.43 

0.90 

0.85 

6.92 

21.16 

17.18 

21.21 

-0.05 

0.90 

0.85 

5.00 

18.42 

16.64 

19.06 

0.64 

0.90 

0.85 

2.00 

13.50 

16.49 

14.10 

0.60 

0,90 

0.85 

1.00 

11.31 

17.35 

11.22 

-0 .09 
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Table  4  (continued) 

Pl/p  2 


/"2 

C1/C2 

.  80 

0.95 

10.0 

.  30 

0,95 

6.92 

.  30 

0.95 

5.00 

1.80 

0.95 

2 . 00 

).30 

0.95 

1.00 

)  .80 

0,90 

10.0 

).ao 

0 ,90 

6.92 

].8a 

0.90 

5 . 00 

0.80 

Q.90 

2.00 

3  .  BO 

0,90 

1.00 

0.30 

0,85 

10.0 

0 . 80 

0 .35 

6.92 

0.30 

0.85 

5.00 

0.80 

0.35 

2.00 

0 . 80 

0,85 

1.00 

0.70 

0.95 

10.0 

0.70 

0,95 

6.92 

0.70 

0,95 

5.00 

0.70 

0.95 

2,00 

0.70 

0.95 

1.00 

0.70 

0.90 

10.0 

0.70 

0.90 

6 ,92 

0.70 

0 .90 

5.00 

0.70 

0,90 

2.00 

0.70 

0 .90 

1.00 

0.70 

0.85 

10.0 

0.70 

0,85 

6,92 

0.70 

0.85 

5,00 

0 .70 

0 . 85 

2.00 

0.70 

0 .85 

1.00 

16.33 

14 . 20 

12.79 

9 . 82 
7 .58 
20.46 
19.14 

16.33 

12.02 

9.82 
25.03 
22.44 
19.14 

14.20 

11.31 
17.07 

14.20 

13 . 50 

9.82 
8.36 
21.16 

19.80 
17.07 

12.79 

10.59 

25.64 
23.12 

19.80 
14.95 
12.02 


17. 5U 

17.21 
17.25 
17.91 
17 . 36 
16.76 

17.70 

16 .80 

16.72 
17.16 
17.47 
17.68 

16.73 

16.83 

16.84 
17.72 

16.64 

17.60 

17.31 

18.51 
16.75 

17.70 
16.98 
17.20 
17.88 
17.30 

17.61 
16.78 
17.13 
17.30 


®DC 
16.03 
14.20 
12.76 
9 . 44 
7  ,  52 
21.02 

18.62 
16 . 54 

12 .33 
9 . 86 
24.67 
21.35 

19 . 64 
14.53 
11.56 
16 . 59 

14.70 

13.21 

9.77 

7.78 
21.75 
19.27 

17 . 31 

12.81 

10.20 

25.52 
22.61 

20.32 
15.03 
11.96 


-0.30 
0 .00 
-0.03 
-0  .  33 
-0.06 
0 .00 
-0 . 52 
0.21 
0 .36 
0 .04 
-0.36 
-0  .  59 
0.50 
0.33 
0.25 
-0.48 
0.50 
-0.29 
0.05 
-0.58 
0.59 
-0  .53 
0.24 
0 .02 
-0 .39 
-0.12 
-0.51 
0.52 
0.08 
0 .06 


Table  4  (continued) 


1/P2 

C1/C2 

®DW 

K 

®DC 

0.50 

0.95 

10.0 

17.74 

16.91 

18.07 

0.33 

0.50 

0.95 

6.92 

15.64 

16 . 82 

16.01 

0 .37 

0.50 

0.95 

5.00 

14.95 

17.90 

14.38 

-0 . 57 

0.50 

0.95 

2.00 

9.82 

15.89 

10 . 64 

0.82 

0.50 

0 .95 

1.00 

9.09 

18.48 

8.47 

-0.62 

0 .50 

0.90 

10.0 

22.44 

16.31 

23.69 

1.25 

0.50 

0.90 

6.92 

21.80 

17 .89 

20.99 

-0 . 31 

0.50 

0.90 

5.00 

18.42 

16.82 

18 . 86 

0 .44 

0.50 

0.90 

2.00 

14.20 

17.53 

13.95 

-0.25 

0.50 

0.90 

1.00 

11.31 

17.54 

11.10 

-0.21 

0.50 

0.85 

10 . 0 

27.20 

16.85 

27  SO 

0 .50 

0.50 

0.35 

6.92 

25.03 

17.50 

24.62 

-0.41 

0.50 

0.35 

5.00 

21.80 

16.96 

22.13 

0.33 

0.50 

0.35 

2.00 

17.07 

17.96 

16.37 

-0.70 

0.50 

0.85 

1.00 

13.50 

17 .34 

13.03 

-0.42 

Average 

K  =  17. 

,22 

Standard 

1  (N-1) 

deviation 

=  0 . 

,  51 

(N-1) 

variance 

=  0  , 

,25 
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APPENDIX  D 


Case  1 

„  ,,  Measured  Measured 


^s  max 

a 

■^DW 

BW 

w 

D 

BW 

2.12  1 

1.60 

24.5 

10.33 

1 

1.10 

22.0 

11.34 

1 

0.33 

19.0 

12.88 

1 

0 .66 

16 . 5 

15.27 

13.5 

15 

1 

0.55 

14.5 

18.38 

15 . 5 

15 

4 

1.60 

23.0 

6.28 

4 

1.10 

22.1 

8.19 

4 

0.83 

19.0 

12.20 

4 

0.66 

16.0 

13.88 

13 . 5 

15 

4 

0  .55 

15.0 

14.71 

15.5 

15 

10  1 

1.60 

17.5 

19 .26 

1 

1.10 

16.0 

20.29 

1 

0.83 

16.5 

20.90 

1 

0.66 

13.5 

22.41 

13 . 5 

15 

1 

0.55 

12 . 5 

23.19 

15.5 

15 

4 

1.60 

16.0 

13.85 

4 

1.10 

13.5 

14.96 

4 

0.83 

13.0 

15.36 

4 

0.66 

12.0 

15.67 

13.5 

15 

4 

0.55 

12.0 

16.33 

15.5 

15 

100  •  1 

1.60 

16.5 

22.68 

1 

1.10 

15.5 

23.27 

1 

0.83 

14.0 

23.91 

1 

0.66 

13.0 

24.81 

13 . 5 

15 

1 

0.  55 

11.5 

26.00 

15  .  5 

15 

4 

1.60 

12.0 

15.18 

4 

1 . 10 

12.0 

14.93 

4 

0.83 

11.5 

15.25 

4 

0  .  66 

11.0 

15.94 

13  .  5 

15 

4 

0.55 

10.5 

16.95 

15.5 

15 

C1/C2  = 

0.8854 

P1/P2 

=  0.8415 

6  =  1.52 
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Case  2 


max 

a 

BW 

w 

Measured 

Measured 

BW 

1 

1.21 

21.5 

13.31 

1 

0.83 

17 . 5 

20.60 

1 

0.63 

14.5 

24.19 

1 

0.50 

12.5 

26.62 

13 . 3 

15.5 

1 

0.42 

10.0 

27.41 

14.0 

14.38 

4 

1.21 

22.5 

15.04 

4 

0.33 

16.0 

19.24 

4 

0.63 

14.0 

19.32 

4 

0.50 

13.0 

20.34 

13.3 

15.5 

4 

0.42 

11.5 

22.75 

14.0 

14.38 

1 

1.21 

18.0 

24.19 

1 

0  .83 

14.5 

25.13 

1 

0 . 63 

14.5 

23.69 

1 

0.50 

12.0 

29.56 

13.3 

15.5 

1 

0 .42 

9.5 

30.18 

14.0 

14.38 

4 

1.21 

16.0 

16.07 

4 

0.83 

14.0 

17.97 

4 

0.63 

14.0 

20.43 

4 

0.50 

12.5 

22.70 

13.3 

15.5 

4 

0.42 

11.5 

24.96 

14.0 

14.38 

1 

1.21 

13.5 

26.64 

1 

0.83 

16.0 

27.72 

1 

0.63 

14.0 

29.29 

1 

0.50 

11.5 

30.89 

13.3 

15.5 

1 

0 . 42 

9.5 

31.55 

14.0 

14.38 

4 

1,21 

14.0 

17.65 

4 

0 .83 

13.5 

17.94 

4 

0.63 

12.5 

19.37 

4 

0  .  50 

11.5 

21.69 

13.3 

15.5 

4 

0 .42 

11.0 

23.80 

14.0 

14.38 

=  0.8854 

P  / D  =  0.3415 

B  =  2.51° 

2 

1  2 
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Case  3 


max 

a 

^DW 

BW 

w 

Measured 

Measured 

3W 

1 

0.58 

13.5 

30.13 

1 

0.40 

8.0 

32.50 

1 

0.30 

7.0 

33.47 

0.24 

7.0 

35.75 

16.4 

25.5 

1 

0.20 

7.5 

38 . 56 

20 

21 

4 

0.58 

14,0 

24.62 

4 

0 . 40 

10.5 

30.65 

4 

0.30 

10,5 

35,06 

4 

0.24 

10.5 

41.51 

15.4 

25.5 

4 

0.20 

11.0 

20 

2 1 

1 

0.  53 

14.5 

24.60 

1 

0.40 

12.5 

30,47 

1 

0  .  30 

11.5 

35.89 

1 

0.24 

11.5 

41.62 

16.4 

25  .  5 

1 

0.20 

12.0 

20 

21 

4 

0.58 

14.5 

33.94 

4 

0.40 

9.0 

35.10 

4 

0  .  30 

7.5 

36.04 

4 

0.24 

7.5 

38.04 

16.4 

25.5 

4 

0,20 

7.5 

40.78 

21 

21 

1 

■X. 

0.58 

14,0 

35,19 

1 

0.40 

8.5 

35.83 

1 

0.30 

7.0 

36.47 

]_ 

0.24 

7.0 

38.37 

16.4 

25.5 

1 

0.20 

7.5 

40.82 

20 

21 

4 

0.58 

13.5 

25.14 

4 

0.40 

11.0 

30.32 

4 

0.30 

10.0 

35.38 

4 

0.  24 

10.0 

40.97 

16 . 4 

25.5 

4 

0.20 

11.0 

35.02 

20 

0  1 

C  X 

C^/C^  =  0.8854  °i'°2  ■  3  =  4.41° 


59 


Case  4 


%ax 

0 

9 

DW 

BW 

w 

Measured 

od 

Measured 

BW 

1 

1.31 

18 . 5 

15.88 

1 

0.90 

14.5 

17.62 

1 

0 . 68 

13.0 

20.75 

1 

0  .  54 

10.0 

22.84 

12 .67 

12.88 

1 

0.45 

8.0 

24,58 

14.7 

13 . 5 

4 

1.31 

17.5 

9.65 

4 

0.90 

15.5 

13 .38 

4 

0.68 

14.0 

17.32 

4 

0.54 

12.0 

21.72 

12 . 57 

12.88 

4 

0.45 

11.0 

25.21 

14.7 

13 . 5 

1 

1.21 

17.0 

21.32 

0.90 

14.0 

22.55 

1 

0.63 

12.0 

24.35 

1 

0  .  54 

10.0 

26.14 

12,67 

12,38 

1 

0.45 

3.0 

26.62 

14.7 

13.5 

4 

1.31 

14.5 

13.80 

4 

0 .90 

13.0 

15.71 

4 

0.68 

12.0 

17.96 

4 

0  .  54 

11.5 

20.56 

12.67 

12,38 

4 

0.45 

10.0 

22.32 

14.7 

13 . 5 

1 

1.31 

16.0 

22.58 

1 

0 .90 

13.5 

23.71 

1 

0.68 

11.5 

25.13 

1 

0.54 

9.5 

26.56 

12 . 67 

12.88 

1 

0.45 

7.5 

26.76 

14.7 

13 . 5 

4 

1.31 

12.5 

15.05 

4 

0 .90 

11.5 

15.87 

4 

0 .68 

11.0 

17.27 

4 

0  .  54 

10.0 

19,32 

12  .67 

12.38 

4 

0.45 

9.0 

21.15 

14.7 

13  .  5 

C^/C2  =  0.9214  '  0-9373  8  =  2.66° 
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CONVERT  ANGLES  IN  DEGREES  TO  RADIANS 
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FROM  THE  BOTTOM  DOMAIN  CALCULATIONS 
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FUNCTION  REFL  (ANGLNI 
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RETURN 

END 


SUBROUTINE  AXI SMR ( QSTR, QF IN t aXLAB,QyL4B , IX,  lY) 

THIS  SUBROUTINE  DRAWS  THE  AXES  AND  WRITES  THE  SC 
MARKINGS  EVERY  ONE-FOURTH  ETA.  THE  VERTICAL  SCAI 
IS  FIXED  SO  VARIOUS  RUNS  CAN  BE  OVERLAID. 
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